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Abstract 

Aim:\Ne performed a comprehensive meta-analysis to determine the association between P2X7 -762T/C polymorphism and 
pulmonary tuberculosis susceptibility. 

Methodology: Based on comprehensive searches of the PubMed, SCI, Elsevier, China National Knowledge Infrastructure 
(CNKI) and Wanfang Database, we identified eligible studies about the association between P2X7 -762T/C polymorphism 
and pulmonary tuberculosis risk. Pooled odds ratio (ORs) and 95% confidence intervals (95%Cls) were calculated in random- 
effects model. 

Results: IK total of 2207 tuberculosis cases and 2220 controls in 8 case-control studies were included in this meta-analysis. 
Allele model (C vs. T: p = 0.15; OR = 0.83, 95% CI = 0.65-1.07), homozygous model (CC vs. TT: p = 0.23; OR = 0.73, 95% 
CI = 0.44 to 1 .22), and heterozygous model (CT vs. TT: p = 0.57; OR = 0.92, 95% CI = 0.68 to 1 .24) did not show increased risk 
of developing pulmonary tuberculosis. Similarly, dominant model (CC+CT vs. TT: p = 0.32; OR = 0.84, 95% CI = 0.59 to 1.19) 
and recessive model (CC vs. CT+TT: p = 0.08; OR = 0.77, 95% CI = 0.57 to 1.04) failed to show increased risk of developing 
pulmonary tuberculosis. Subgroup analysis by ethnicity did not detect any significant association between P2X7-762T/C 
polymorphism and pulmonary tuberculosis susceptibility. 

Conclusions: P2X7 -762T/C gene polymorphism is not associated with pulmonary tuberculosis susceptibility. 
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Introduction 

Tuberculosis is caused by the bacillus Mycobacterium tuberculosis 
and remains a major challenge to global public health [1]. 
According to the report of World Health Organization, there were 
an estimated 8.6 million incident case of tuberculosis and 1.3 
million deaths in 2012 [2] . Although one-third of population is 
infected by Mycobacterium tuberculosis^ only 10% of those develop 
clinical disease during their lifetime. Multiple factors contribute to 
the risk of infection and development of tuberculosis including 
environmental factors, host-pathogen interactions and genetic 
factors [3]. Epidemiological studies indicate that the risk of 
developing tuberculosis in human is strongly influenced by genetic 
factors [4]. So far, gene polymorphisms of SLCllAl (formerly 
JVRAMPl), vitamin D receptor, toll-like receptor 2, tumor necrosis 
factor-alpha and monocyte chemoattractant protein- 1 have been 
identified to be associated with susceptibility of tuberculosis [5-9] . 

Human P2X7 gene encoding the P2X7 receptor contains 13 
exons and is localized on chromosome 12q24 [10]. P2X7 
receptors are abundantly expressed in macrophages. Mycobacterium 
tuberculosis is a facultative intracellular pathogen and macrophages 
are the main reservoir for this agent. Activation of P2X7 by 



adenosine triphosphate (ATP) causes an immediate opening of a 
cation selective channel, permitting the influx of Ca^^ and Na^ 
and the efflux of K^. This process results the induction of the 
caspase cascade, apoptosis, and the activation of phospholipase D. 
Phospholipase D promotes phagosome-lysosome fusion and 
causes mycobacterial death [11-13]. 

Previous studies reported that various single-nucleotide poly- 
morphisms (SNPs) in P2X7 gene lead to the loss of receptor 
function, and the most common SNPs are the 1513A/C and - 
762T/C [14]. P2X7 1513 A/C gene polymorphism has been 
reported to be significantly associated with increased tuberculosis 
susceptibility [15]. However, the results from different studies on 
the association between the P2X7 -762T/C (rs2393799) polymor- 
phism and pulmonary tuberculosis susceptibility are conflicting 
[16-20]. Meta-analysis is a useful method for combining the 
results from different studies and providing more reliable 
conclusions [21]. Therefore, we performed a meta-analysis to 
investigate the association between the P2X7-762T/C gene 
polymorphism and the risk of pulmonary tuberculosis. 
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Materials and Methods 

Literature Search Strategy 

We conducted a systematic search on the PubMed, SCI, 
Elsevier, China National Knowledge Infrastructure (CNKI) and 
Wanfang Database (up to December 12, 2013). The following 
terms were used: ("tuberculosis" or : "TB") and ("Pi'X/" or 
"P2iL¥7") and ("polymorphism" or "mutation" or "variation"). 
We also searched the references of the retrieved literatures for 
additional studies. There was no restriction for time period, sample 
size, population, language, or types of reports. If more than one 
article was published using the same case series, only the study 
with largest sample size was included. 

Inclusion and Exclusion Criteria 

Included studies in this meta-analysis met the following criteria: 

(a) a human case-control study on the association between P2X7 - 
762T/C polymorphism and the risks of pulmonary tuberculosis; 

(b) containing available genotype data in cases and controls for 
estimating an odds ratio (OR) and 95% confidence interval (CI); 

(c) genotype distributions of control population were consistent 
with Hardy- Weinberg equilibrium (HWE). The exclusion criteria 
were: (a) reviews, letters, editorial articles and case reports; (b) 
studies on the association between other gene polymorphisms and 
pulmonary tuberculosis risks. 



272 studies were searched from PubMed, 
SCI, Elsevier, China National Knowledge 
Infrastructure (CNKI) and Wanfang 
Database 








256 Studies were excluded 
after title or abstract review 








16 Studies chosen for full text review 








5 Studies were excluded for 
not relevant to -762 T C 
polymorphism 




11 Studies screened for a meta-analysis 








3 Studies were excluded: 1 
was overlapped 1 was not 
consistent with HWE, and 1 
without sufficient usable 
data 




8 studies included m qualitative s\-nthesis 
(meta-analysis) 





Figure 1. Flow diagram of the study selection. 

doi:1 0.1 371 /journal. pone.0096359.g001 



Data Extraction 

Two investigators (Yi and Cheng) extracted the data from all of 
the eligible publications according to the inclusion and exclusion 
criteria. Primary extraction data were reviewed by Zhen, and any 
disagreement was resolved by discussion among the authors. From 
each study, the following information was extracted: first author's 
name, year of publication, study location, sample size, source of 
control, the genotyping method, the number of genotype 
frequencies in cases and controls. 

Statistical Analysis 

The association between P2X7 -762T/C polymorphism and 
pulmonary tuberculosis risks was estimated by calculating pooled 
ORs and 95%CI in the allele model (C vs. T), co-dominant model 
(CC vs. TT and CT vs. TT), dominant model (CG/CT vs. TT), 
and recessive model (CC vs. CT/TT). Significance of the pooled 
ORs was determined by Z test (p<0.05 was considered 
significant). Heterogeneity among included studies was checked 
by chi-square-based Qtest and test [22]. If the data showed no 
heterogeneity (p>0.05, I^<50%), Mantel-Haenszel fixed effects 
model was used, otherwise DerSimonian-Laird random effects 
model was used [23]. Hardy- Weinberg equilibrium (HWE) was 
also tested for included studies. Sensitivity analyses were 
performed by omitting certain studies each time. Publication bias 
was assessed by Begg's adjusted rank correlation test and Egger's 
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Experimental 



Control 



Odds Ratio 



Study or Subgroup 




Total Events 


Total 


Weight 


M-H 


Random. 95% CI 


Bahari 2013 


196 


300 


248 


300 


11.4% 




0.40 [0.27. 0.58) 


Ben-Selma 2011 


89 


336 


79 


300 


11.9% 




1.01 [0.71, 1.43] 
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Minr>-Mr>r#ann 9007 


48 


184 


74 


220 


10.6% 




0 70 fO 45 1 071 


Sambaslvan 2010 


164 


312 


79 


200 


11.8% 




1.70 [1.18. 2.43) 


Singia 2012 


401 


572 


605 


784 


13.6% 




0.69 [0.54, 0.89) 


Songane 2012 


775 


1684 


766 


1688 


15.0% 




1.03 [0.90, 1.18) 


Total (95% CI) 




4414 




4440 


100.0% 




0.83 [0.65,1.07] 


Total events 


2096 




2255 










Heterogeneity: Tau^ = 


0.10: Chi^ 


= 43.87. df 


= 7(P 


< 0.00001); P = 


84% 





Odds Ratio 
M - H . Random. 95% CI 



Test for overall effect: Z = 1 .45 (P = 0.15) 



0.01 0.1 1 10 100 

Favours [experimental) Favours [control) 



Figure 2. Forest plot and ORs with 95% CI of P2X7-762T/C polymorphism and pulmonary tuberculosis risk (C vs. T). The squares and 
horizontal lines correspond to the OR and 95% CI for each study. The area of the squares reflects the weight. The diamond represents the summary 
OR and 95% CI. 

doi:10.1371/journal.pone.0096359.g002 



regression asymmetry test.p<0.05 was considered as statistically 
significant [24]. All statistical analyses were performed using 
RevMan 5.0.23 (Cochrane Library Software, Oxford, UK) and 
STATA 1 1.0 Software (StataCorp LP, College Station, TX, USA). 

Results 

Studies Selection Process and Characteristics 

A total of 272 studies were retrieved through database 
searching. After reading the titles and abstracts, 16 potential 
studies were identified for further investigation. Five studies were 
excluded for not being relevant to -762 T/C polymorphism. One 
study was excluded for data overlapping and one study was 
excluded for not reporting the usable data. One case-control study 
was excluded for being not consistent with HWE. Finally, 8 studies 
[16-18,20,25-28] were included for data extraction (Figure 1). 
Among these studies, 6 studies focused on the association between 
pulmonary tuberculosis risks and P2X7 -762T/C polymorphism 
while other studies enrolled both pulmonary tuberculosis and 
extra-pulmonary tuberculosis patients. All of the 8 included studies 
were non-familial case-control studies. There was a total of 2207 
subjects with pulmonary tuberculosis 1442 subjects were diag- 



nosed by culture, 479 subjects were diagnosed by smear, and the 
methods of diagnosis of the rest subjects were not reported 
(Table 1). Genotype and allele distributions for each case-control 
study are shown in Table 2. 

Quantitative Data Synthesis 

A total of 2207 cases and 2220 controls in 8 case-control studies 
were pooled together for evaluation of the overall association 
between P2X7 -762T/C polymorphism and risk of pulmonary 
tuberculosis. The pooled OR from all studies indicated no 
significant association between P2X7 -762T/C polymorphism 
and pulmonary tuberculosis risk in allelic model (C vs. T: p = 0. 15; 
OR = 0.83, 95% CI = 0.65-1.07) (Figure 2). There was no 
significant association in homozygous model (CC vs. TT: 
p = 0.23; OR = 0.73, 95% CI = 0.44 to 1.22) or heterozygous 
model (CT vs. TT: p = 0.57; OR = 0.92, 95% CI = 0.68 to 1.24) 
(Figure 3, 4). Similarly, dominant model (CC+CT vs. TT: 
p = 0.32; OR = 0.84, 95% CI = 0.59 to 1.19) and recessive model 
(CC vs. CT+TT: p = 0.08; OR = 0.77, 95% CI = 0.57 to 1.04) did 
not demonstrate any altered risk for pulmonary tuberculosis 
(Figure 5, 6). 
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99 
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M-H. Random. 95% CI 



Test for overall effect; Z = 1.21 (P = 0.23) 
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Figure 3. Forest plot and ORs with 95% CI of P2X7 -762T/C polymorphism and pulmonary tuberculosis risk (CC vs. CT-I-TT). 

doi:10.1371/journal.pone.0096359.g003 
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Figure 4. Forest plot and ORs with 95% CI of P2X7 -762JIC polymorphism and pulmonary tuberculosis risk (CC vs. TT). 
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Figure 5. Forest plot and ORs with 95% CI of P2X7 -762T/C polymorphism and pulmonary tuberculosis risk (CT vs. TT). 
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Subgroup analysis by etlinicity did not detect any significant 
association between P2X7-762 T/C and pulmonary tuberculosis 
risk in Asians or Caucasians. A summary of results from all 
comparisons is listed in Table 3. 



Publication Bias 

Begg's funnel plot and Egger's test were used to assess the 
publication bias of included studies. Publication bias was not 
observed in Begg's funnel plot. The shape of the funnel plots was 
symmetrical and the Egger's test did not show any evidence of 
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Figure 6. Forest plot and ORs with 95% CI of P2X7 -762T/C polymorphism and pulmonary tuberculosis risk (CC+CT vs. TT). 

doi:10.1371/journal.pone.0096359.g006 
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Table 3. Meta-analysis of P2X7 gene -I^IJ/C polymorphism and risk of pulmonary TB in each subgroup. 





Category 


C vs. T 


CC vs. TT 


CT vs. TT 


CC+CT vs. TT 


CC vs. CT+TT 


OR (95%CI), 1^ (%) 


OR (95%CI), P (%) 


OR (95%CI), P (%) 


OR (95%CI), P (%) 


OR (95%CI), 1^ (%) 


Ethnicity 


Caucasian 


0.86 (0.60-1.24), 89 


0.79 (0.37-1.68), 87 


0.97 (0.57-1.67), 77 


0.85 (0.45-1.61), 86 


0.81 (0.56-1.19), 79 


African 


0.82 (0.57-1.17), 66 


0.66 (0.31-1.41), 62 


0.82 (0.63-1.08), 0 


0.80 (0.56-1.16), 48 


0.69 (0.37-1.31), 50 


Mixed 


0.70 (0.45-1.07) 


0.52 (0.20-1.34) 


0.71 (0.39-1.30) 


0.66 (0.38-1.16) 


0.60 (0.24-1.49) 


Overall 


0.83 (0.65-1.07), 84 


0.73 (0.44-1.22), 81 


0.92 (0.68-1.24), 67 


0.84 (0.59-1.19), 79 


0.77 (0.57-1.04), 69 
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publication bias. These data indicate that there is no significant 
publication bias in this meta-analysis. 

Test of Heterogeneity 

Heterogeneity among studies in each model is shown in Table 3. 
The showed a stable variation under all comparisons (C vs. T: 
1^ = 84%; CC vs. TT: 1^ = 81%; CT vs. TT: 1^ = 67%; CC vs. 
CT+TT: 1^ = 79%; CC+CT vs. TT: 1^ = 69%). Hence, the 
random efTect model was applied to calculate the pooled ORs 
and 95% CI. In the subgroup analyses, p value for heterogeneity 
was not significant under heterozygous, dominant and recessive 
models of Asian population. However, there were heterogeneity 
under all comparisons in Caucasian population (Table 3). 

Sensitivity Analysis 

Sensitivity analysis was performed to assess the influence of each 
individual study on the pooled OR by deleting one single study 
each time. The results showed that no individual study affected the 
pooled OR significantly, indicating the stability of this meta- 
analysis. 

Discussion 

In the present meta-analysis, we show that P2X7 gene -762T/ C 
polymorphism is not associated with the risk of pulmonary 
tuberculosis. Subgroup analysis by ethnicity did not detect any 
significant association between P2X7-1^2 T/C and pulmonary 
tuberculosis risk in Asians or Caucasians. 

Our meta-analysis included a total of 8 case-control studies with 
2207 cases and 2220 controls, which is substantially larger than a 
previous meta-analysis on P2X7 gene -762T/C polymorphism and 
pulmonary tuberculosis risk [15]. We performed the comparisons 
of all genetic models in this meta-analysis, while only allelic 
frequencies were compared in the previous meta-analysis [15]. 
The conclusion of our meta-analysis that P2X7 gene -762T/C 
polymorphism is not associated with pulmonary tuberculosis risk is 
consistent with the previous meta-analysis [15]. 

The polymorphisms in the promoter or coding region of P2X7 
gene may alter its expression or function. However, the -762T/C 
polymorphism of P2X7 does not appear to be related with altered 



receptor expression [16-18]. Li and colleagues found that patients 
with active tuberculosis exhibited reduced P2X7 receptor expres- 
sion levels on mononuclear cells which were recovered during 
chemotherapy. This suggested that the expression of P2X7 
receptor may be regulated by host- or pathogen-generated 
inhibitory factors [16,25,26]. 

Heterogeneity is a potential problem that might affect the 
interpretation of the results. Therefore, the control group should 
conform to HWE to ensure similarity of the population genetic 
background and the reliability of association analyses. One of the 
potentially eligible studies did not conform to HWE and was 
excluded [19]. The heterogeneity could also be caused by ethnic 
variations, environmental factors and methodological factors in 
study design. The leave-one-out sensitivity analysis did not affect 
the results of this pooled analysis, indicating that our results were 
reliable. There was no publication bias for the association between 
P2X7-762 T/C polymorphism and pulmonary tuberculosis risk. 

The present meta-analysis had several limitations. First, we only 
included the studies published in English and Chinese. Some 
related studies published in other languages might be missed in our 
meta-analysis. Second, although the funnel plot and Egger's test 
did not show any publication bias, the influence of bias in the 
present analysis could not be completely excluded. A more precise 
analysis could be performed if more detailed data such as age, sex 
and exposure of each subject was available. 

In conclusion, this meta-analysis indicates that P2X7 -762T/C 
polymorphism is not associated with pulmonary tuberculosis risk. 
Further studies with more stringent design and larger sample size 
are required to validate this conclusion. 
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